Anisotropic thermally activated diffusion in percolation systems 
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We present a study of static and frequency- dependent diffusion with anisotropic thermally acti- 
vated transition rates in a two-dimensional bond percolation system. The approach accounts for 
temperature effects on diffusion coefficients in disordered anisotropic systems. Static diffusion shows 
an Arrhenius behavior for low temperatures with an activation energy given by the highest energy 
barrier of the system. From the frequency-dependent diffusion coefficients we calculate a character- 
istic frequency luc ~ f/ic, related to the time tc needed to overcome a characteristic barrier. We 
find that lJc follows an Arrhenius behavior with different activation energies in each direction. 

PACS numbers: 05.40.-a 05.60.Cd 66.30.-h 



The study of diffusion on disordered media is an im- 
portant problem, in view of its relevance in a wide variety 
of natural and industrial processes 0, 0, 111 • In the last 
years the anisotropic generalization of diffusion has at- 
tracted much attention (J, |5j J,, ^ ^] , justified by the 
diversity of systems in which diffusion takes place. A 
few examples of anisotropic systems are porous reservoir 
rocks 3, 8, 10], epoxy-graphite disk composites and 
layered semiconducting compounds [l2j |. 

Static (long-time) and frequency-dependent conductiv- 
ity on isotropic disordered media has been extensively 
studied with both analytical and numerical methods. 
One of the most widely used models for disorder me- 
dia is the percolation model 0, 0, 0| , due to its sim- 
plicity and interesting properties (characteristic percola- 
tion threshold, fractality of the sample-spanning cluster, 
etc). More recently, diffusion has also been studied in 
anisotropic bond percolation systems 0,'5', However, 
the interplay between temperature and disorder in these 
systems has not been studied yet. The natural manner 
to include temperature effects in diffusion problems is 
through thermally activated processes, i.e., associating 
energy barriers to the transition rates. Interesting re- 
sults are found at low temperatures in isotropic systems, 
where the competing effects of temperature, energy bar- 
riers and topology become important [la, UgL . Static 
diffusion is then described by an Arrhenius law with a 
characteristic energy which depends on the percolation 
threshold of the lattice 0, , and frequency-dependent 
diffusion becomes universal if properly scaled [l3- The 
scaled units include information about different parame- 
ters: temperature, characteristic percolation energy, and 
a characteristic frequency that marks the onset of static 
diffusion. 

In this paper we shall focus on the description of a 
model for anisotropic diffusion processes, both in the 
static and frequency-dependent regimes. In order to em- 
phasize the role of temperature and its relevance for diffu- 
sion in a system with energy disorder, we use here a two- 
dimensional isotropic bond percolation system with dif- 
ferent energy-dependent intrinsic transition rates in each 



direction. We define the transition rates as Wi(2) in the 
1(2) direction and choose them according to r(wi) and 
^1{'W2), the probability distribution functions (PDFs), 
which for this anisotropic model are given by 



r(u;i) : 



pS{wi ^ 

-- p6{'W2 



(1- 
(1 



-p)6iwi), 
-p)5{w2). 



(1) 



It means that the transition rate wi takes the value Wi 
with occupation probability p, and zero otherwise, and 
analogously for the 2 direction. 

In order to account for the temperature dependence 
of diffusion coefhcients, we propose intrinsic transition 
rates characterized by a thermally activated process with 
different activation energies in each direction. Therefore 
we define an anisotropic thermally activated process, in 
which the transition rates in each direction take the form 



^1(2) 



7oexp 



( gi(2)\ 
V kT J 



(2) 



where 70 is the constant jump rate and k is the Boltz- 
mann constant. The (1 — p) fraction of non-conductor 
components in Eq. Q now represents the existence of 
infinite energy barriers. 

We set ei > £2 and define an anisotropic parameter 
a = ei/e2, and a mean energy e = (ei -I- e2)/2. This 
mean energy was kept constant in the present work, fn 
terms of these parameters we may write ei = 2ae/{a + 1) 
and £2 — 2£/(q! + 1). Then, the relevant parameters of 
the problem become the occupation probability p, the 
temperature T, and the anisotropic parameter a. In the 
following, energies and temperatures are measured in ar- 
bitrary units (with fc = 1). 

The proposed model is studied both analytically, by us- 
ing an anisotropic extension to the effective medium ap- 
proximation (EM A), and numerically by means of stan- 
dard Monte Carlo (MC) simulations. We shall briefly 
describe both methods. 

The reader is referred to jSj) for a complete descrip- 
tion of anisotropic EMA. Here we only summarize the 
key results. The EMA consists in averaging the effects 



2 



of disorder by defining an eff'ective medium with eff'ective 
transition rates, which depend on the Laplace variable u. 
These effective transition rates are self-consistently deter- 
mined by the requirement that the difference between the 
propagator of the impurity and homogeneous problems 
should average to zero. Thus, in anisotropic problems, 
two effective transition rates, w\{u) and ^^(u) (one for 
each direction), are introduced. These effective transi- 
tion rates are determined by two self-consistent condi- 
tions 



wi — Wi 



W2 — W2 



= 0, 



= 0. (3) 



Here, G^'^^ and are the non-pcrturbed anisotropic 
Green functions related to the probabilities of moving 
from the origin to one of its nearest neighbors in the 
1(2) direction and the return probability, respectively. 
The impure bond connects two nearest neighbor sites of 
the lattice whose transition rates are equal to wi if the 
impure bond lies in the 1 direction and W2 if the impure 
bond is in the other direction. 

Following the linear response theory the general- 
ized frequency-dependent complex diffusion coefficients 
D[u)) in the anisotropic EMA context are given by 
Di(2)(.'^) — ^1(2) ('^ ~ *'^)' where a is the lattice con- 
stant. In the following we will take a — 1. 

In the static case [uj — 0) Eqs. Q become 



Dl - 70 exp 



kTJ 



arctan 



-|-P70 exp 



70 exp 



V kTJ 



arctan 




-fP7oexp — —\ 



kTJ 



= 0, 



= 0, (4) 



where the transition rates are defined in Eq. (0), D\ = 
Di{uj = 0) = wf(0) and = D2{cj = 0) = w^{0) 
represent zero-frequency (static) diffusion coefficients in 
each direction. 

Monte Carlo data are obtained by performing clas- 
sical random walk simulation in square lattices with 
300^ sites. Normal diffusion regime (long-time limit) is 
reached between 10'^ and 10® steps, and the mean square 
displacement is averaged over 2000 and 10000 different 
lattices, depending on data fluctuations. Static diffu- 
sion coefficients are obtained via the Einstein's relations, 

(-^1(2)) = 2-^1(2)*- '^^^ shown earlier standard 
MC simulations fail to describe particle diffusion for very 
low temperatures. In this limit the method is inefficient 
and other simulation methods are necessary. Here the 
standard MC method was used and simulations were per- 
formed for temperatures above T = 0.05. 
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FIG. 1: Dependence of the anisotropic static diffusion coef- 
ficients on the occupation probability p for T = 0.125. The 
characteristic value p* is also indicated. Energy and tempera- 
ture are given in arbitrary units with the Boltzmann constant 
= 1. 



We present now the results for the long-time diffu- 
sion properties. Firstly, in the isotropic case a = 1, 
the transition rates are given by 70 exp (—j^) and, as 
expected, only one diffusion coefficient is obtained, viz.. 



^2 



70 exp {-^) (2p - 1). Secondly, 
for all values of a, in the T — > 00 limit (kT ^ e) 
the model again reduces to the isotropic bond perco- 
lation problem with transition rates 70 [Eq. In 
this case, the isotropic diffusion coefficient is given by 
-^iso — 7o(2p — 1). Finally, for values of a 1 and in 
the low temperature limit and D2 were calculated by 
numerically solving the set of Eqs. Q), and results are 
compared with MC simulations. 

Let us consider the variation of and D2 with the oc- 
cupation probability p, for fixed values of T in anisotropic 
conditions. The calculation was performed for e = 0.5 
and 7o — 0.25j^, where v is a, characteristic jump fre- 
quency. Figure ^ shows -D5'(2) ^ function of p for 
T — 0.125. Symbols represent MC simulations and lines 
correspond to the EMA numerical solution of Eqs. I@J. 
The isotropic a — 1 case is included for comparison. For 
a > 1, we obtain < D2, which is in agreement with 
the fact that lower energy barriers imply higher diffu- 
sion coefficients. For p < pc = 0.5, the isotropic bond 
percolation threshold [13, we find that Di — D2 — 
as expected. In the low p region the two anisotropic 
diffusion coefficients are lower than the isotropic one 
(p < p* — 0.68). In this region the conducting cluster 
is poorly connected and the diffusion in the low energy 
direction is highly affected by the existence and height of 
high energy barriers. As p is increased, more energy bar- 
riers appear and the diffusion in the low energy barrier 
direction is less sensitive to high energy barriers. This 
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FIG. 2: Arrhenius plot for the anisotropic static diffusion 
coefficients at low temperatures for p — 0.8. Energy and 
temperature are given in arbitrary units with the Boltzmann 
constant k — 1. 



kind of behavior is also found for other values of T, with 
increasing values of p* as T decreases. 

Next we turn to the behavior of and D2 at low tem- 
peratures for fixed values of p. We do not present here 
results for p = pc (or near pc), the percolation threshold, 
where the substrate has fractal properties and diffusion 
becomes anomalous 0, • In Fig- 121 we show the tem- 
perature dependence of the static diffusion coefficients in 
an Arrhenius plot for p = 0.8 and a = 1 and 2. From the 
EMA curves (lines) the slope of the linear region is nu- 
merically calculated. For isotropic media, the EMA pre- 
dicts a slope equal to —0.5 = — e. For a = 2, we find that 
Di (T) and D2 (T) reach an asymptotic Arrhenius behav- 
ior with the same activation energy ei, corresponding to 
the highest energy of this anisotropic model. This is not 
a striking result considering that long-time diffusion at 
low temperatures is strongly governed by the highest en- 
ergy barrier, as the particle spends a lot of time trying 
to overcome them while performing long-distance trajec- 
tories. The low temperature behavior of static diffusion 
coefficients depends only on the highest characteristic en- 
ergy value of this anisotropic model. 

The frequency behavior of diffusion coefficients was 
studied by numerically solving the set of coupled 
Eqs. JSJ, using the corresponding expressions for the 
anisotropic Green functions 5] , and the temperature de- 
pendence introduced in Eq. |(5J. We present here re- 
sults for the real part of the frequency-dependent com- 
plex diffusion coefficients. In Fig. |31we show the scaled 
values iDi(2)(w)/-D^2) fo^' P — ^-8, e — 0.5, and temper- 
ature values between 0.1 and 0.04. Three characteristic 
regimes can be distinguished: (?) the long-time limit for 
w — > 0, (ii) a power law behavior characteristic of inter- 
mediate frequencies, a ^ lo'^, with s < 1 and, (iii) the 
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FIG. 3: Scaled frequency-dependent diffusion coefficients for 
isotropic (a) and anisotropic a — 2 [(b) and (c)] cases. Each 
curve is labeled with its corresponding temperature (in arbi- 
trary units). 



high-frequency regime, lu — > 00, where the diffusion coef- 
ficients approach to a constant value. This kind of dielec- 
tric response has been observed before in a broad class of 
ionic and electronic systems, in isotropic and anisotropic 
media [HiSIIll. 

In the isotropic case we find that while decreasing 
the temperature the scaled value Diso{i~^)/ D^^o reaches 
a temperature-independent value for w ^ cx) [Fig. |3| 
(a)]. In the anisotropic case, while Di{uo)/D\ satu- 
rates for w — > 00 to a temperature-independent value 
at low temperatures [FigOKc)], D2{uj)/D2 saturates to a 
temperature-dependent value that follows an Arrhenius 
law with an activation energy equal to ei — £2, even at 
low temperatures [FigOKb)]. This lo ^ 00 behavior can 
be interpreted from the evolution of the particle diffu- 
sion for i — > 0. For a = 1 we find that Diso{ijJ — > 00) = 
P7o exp (— . This means that the first step of the par- 
ticle is given by the probability p of finding a given bond 
times the intrinsic transition rate of that bond. From this 
considerations and using the isotropic result for D^^^{T) 
we find that Aso(w 00) /D%^ = p/{2p - 1) = 4/3 
for p = 0.8, as shown in Fig. |21 (a). In anisotropic 
conditions, we expect the < — > (tj — > 00) evolution 
to be given by p7oexp(— ^) in the 1 direction and 
by p7oexp(— -ly) in the 2 direction. We showed pre- 
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FIG. 4: Temperature dependence of the characteristic fre- 
quencies u)^n2)- The isotropic case uoc,iso (full squared) and 
the anisotropic q = 2 cases Uci (open circles) and luc2 (open 
diamond) are shown. Lines are guides to the eye and are la- 
beled by its corresponding activation energy. The inset shows 
the determination of the value luc from a typical curve for the 
frequency-dependent diffusion coefficient. The value m rep- 
resent the slope of the tangent curve through the inflection 
point {\og-^Qijj* /v,D* /a^v) (see text). 

viously that D\ ^ exp(— -1^), and Z?" ^ '^^P {~It) ^ 
thus Di{ui oo)/Di is expected to be constant at low 
temperatures, and D2{uj — + oo)/Z?§ ^ exp (^^j^), which 
account for the calculated EMA results for cj ^ oo in Fig. 
13 (b) and (c). 

For intermediate frequencies a characteristic frequency 
ujc ^ 1/tc can be defined, related to the time tc needed 
to overcome a characteristic energy barrier. The inset 
of Fig. ^shows a typical frequency-dependent D{ui) plot 
(as those in Fig. O. The ujc parameter is calculated from 
our anisotropic EMA data as follows. Let the inflection 



point be (log^o and m the slope of the 

tangent curve passing through the inflection point. The 
^ci(2) parameter is given by the intersection of the slope 
m through the inflection point and the line corresponding 
to the long-time diffusion coefficient D^^) (Fig- Bl inset). 
We include the subscript 1(2) as lUc may be different for 
each direction. In Fig. ^ we present the temperature 
dependence of the characteristic frequencies Wci(2) in an 
Arrhenius plot for a — 1 and 2. All cases reach a linear 
behavior for low temperature values, but different slopes 
are found in each case, corresponding to the characteris- 
tic energy in each direction. This indicates, as expected, 
that different times ici(2) exp ( — ^r-) are needed in 
each direction to overcome its characteristic energy bar- 
riers. 

In summary, in the present paper we described long- 
time and frequency-dependent diffusion in anisotropic 
thermally activated processes in a two-dimensional bond 
percolation lattice. The present model has different ac- 
tivation energies in each direction of the square lattice. 
Long-time diffusion coefficients follow an Arrhenius law 
at low temperatures with the highest energy barrier be- 
ing the activation energy for diffusion in both directions 
of the lattice. From the frequency-dependent diffusion 
coefficients, we define a characteristic frequency uJci{2)- 
We remark that Uc is not marking the onset of long-time 
diffusion, as in isotropic problems with continuous PDFs 
[TtI . but it is just the frequency associated to overcome 
the characteristic energy of each direction. The exis- 
tence of two different activation energies for a;ci(2) and 
only one for D^^^ is a consequence of the strength of the 
percolation-type disorder, in contrast with results from 
isotropic continuous PDFs, where the characteristic fre- 
quency and the static diffusion coefficient follow Arrhe- 
nius laws with the same activation energy 0, 0, 0| . 
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